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It was proposed before that single bubble sonoluminescence (SBSL) may be caused by strong 
electric fields occurring in water near the surface of collapsing gas bubbles because of the flexoelectric 
effect involving polarization resulting from a gradient of pressure. Here we show that these fields can 
indeed provoke dynamic electric breakdown in a micron-size region near the bubble and consider the 
scenario of the SBSL. The scenario is: (i) at the last stage of incomplete collapse of the bubble, the 
gradient of pressure in water near the bubble surface has such a value and a sign that the electric field 
arising from the flexoelectric effect exceeds the threshold field of the dynamic electrical breakdown 
of water and is directed to the bubble center; (ii) mobile electrons are generated because of thermal 
ionization of water molecules near the bubble surface; (iii) these electrons are accelerated in "cold" 
water by the strong electric fields; (iv) these hot electrons transfer noble gas atoms dissolved in water 
to high-energy excited states and optical transitions between these states produce SBSL UV flashes 
in the trasparency window of water; (v) the breakdown can be repeated several times and the power 
and duration of the UV flash are determined by the multiplicity of the breakdowns. The SBSL 
spectrum is found to resemble a black-body spectrum where temperature is given by the effective 
temperature of the hot electrons. The pulse energy and some other characteristics of the SBSL are 
found to be in agreement with the experimental data when realistic estimations are made. 
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I. INTRODUCTION 



Sonoluminescence refers to the phenomenon of light emission during acoustic radiation of a liquid and is associated 
with cavitation bubbles present in the liquid. The most controllable and promising experimental data were obtained 
for single bubble sonoluminescence (SBSL): picosecond UV radiation of a single bubble pulsating in the field of the 
sound wave . 

Though the case of a single pulsating bubble is, of course, much simpler than cavitation in general, it turns out that 
the SBSL is a very complex phenomenon which still remains not completely understood. There is a vast literature 
devoted to SBSL that was recently reviewed in Ref. Interest in this phenomenon is stimulated, on the one hand, 
by the fact that in relatively simple and controllable experiments extraordinary conditions (ultra-high pressures, 
temperatures, and ultra-short light flashes) are realized at the final stage of the bubble collapse; on the other hand, 
it looks promising to use the SBSL to construct a source of ultraviolet ultra-short flashes that is much cheaper than 
lasers. 

The highly involved hydrodynamics of bubble collapse has been addressed in many papers that are reviewed in Ref. 
. One of the most important problems to explain here is the very existence of stable pulsation of bubbles. Recently 
a theory has been developed which explains this regime in terms of the dissolved gas diffusion and chemical 
reactions in the gas within the bubble. It has been shown that an accumulation of noble gas in the bubble takes place 
during pulsation of the bubble and that in a stable situation the gas in the bubble consists, almost entirely, of noble 
gas and, of course, of water vapor. These theoretical conclusions were supported by the experimental data of Ref. [||. 

Another key problem in SBSL is the mechanism of the light emission. Many mechanisms have been proposed, 
criticized, and reviewed to explain UV flashes radiated by the collapsing bubble ||. The most popular of the 
suggested mechanisms are the adiabatic heating of the bubble gas > shock wave-Bremsstrahlung model |p|,puj|TT2| , 

and, recently proposed proton-tunneling radiation as a result of a phase transition in water at ultra-high pressure 
jL3|| . A general feature of these mechanisms is that extraordinary conditions are needed which can only be realized, 
if at all, at a very small bubble radius when the density of the bubble gas is close to the water density and both the 
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hydrodynamics of the bubble and the properties of the bubble content and the neighboring liquid are not known from 
experiment but inferred from numerical simulations. 

In this paper we show that the main characteristic features of SBSL can be explained even without making any 
assumptions about the extraordinary conditions. We are far from saying that these extraordinary conditions are 
not present in experiments. What we mean is that there is another mechanism of the SBSL which occurs in water 
near the bubble surface but not in the bubble gas as has been assumed in the most popular models of SBSL ||. 
Extraordinary conditions are not necessary for the operation of this mechanism even though it may well be the main 
cause of UV radiation. The mechanism under consideration is based on the idea put forward in Ref. |l4| that SBSL 
occurs because of electric breakdown in strong electric fields arising near the bubble surface as a result of flexoelectric 
effect, that is the effect of polarization of water because of gradients of pressure fll6f . Here we present a scenario 
and estimations that show that within the hypothesis the main features of SBSL can be explained using relatively 
moderate parameters, e.g., temperatures (5-^10)10 3 K for the bubble gas and the natural (expected) value of the 
flexoelectric coefficient of water jl4| . 

Let us mention that Ref. [Q has left many questions to answer. The origin of the optical radiation of the bubble 
remaines unexplained. The problem is that in the visible region pure water has very low levels of absorbtion and 
radiation due to interband optical transition llTfl. Moreover the breakdown scenario is far from being clear. The 
mechanism of the breakdown in water (see, e.g., ||18[|) involves "lucky electrons" whose acceleration in the electric field 
leads to development of an avalanche. At ambient temperature the concentration of free electrons in water is quite 
negligible (the band gap is about 6.5 eV [jl7|,[l8)) and it is impossible to find a "lucky electron" in a small volume near 
the bubble surface, i.e. in the region of high electric field, during the short time of the existence of this field. The 
spectacular synchronization of the emission pulses [§,[l9| that in the case of breakdown appears, at least at first sight, 
to be hardly compatible with the fact that we are dealing with a probabilistic situation. In addition, the reference 
to the Penning effect to explain the role of the noble gases is not convincing because the ionization energy of the 
metastable state is more than the band gap of water as distinct from the case of breakdown in gases where the Penning 
effect is pronounced. 

In this paper we consider in more detail possible processes associated with the electric field arising from the 
flexoelectric effect near a bubble exhibiting SBSL. A common difficulty with theories of SBSL is that little is known 
about the parameters of gas in the bubble when the bubble radius is near its minimum value. Consensus exists about 
one point: it is far from being a gas, because the lowest limit of the radius is governed by the van der Waals repulsion 
and the minimum volume of the bubble is close to the van der Waals hard core volume. The equations of state used 
for these conditions are not reliable at present. This is a challenging and fundamental problem but it is beyond the 
scope of this paper. What is now possible is to make order-of-magnitude estimations. That is why we shall first 
explain what values of relevant parameters are necessary for our scenario to be operative and then discuss whether 
our assumptions about these parameters are realistic. We include the value of the flexoelectric coefficient of water 
among these parameters as well. This coefficient has not been measured, unfortunately, and we use its estimated 
"natural" value but, from the other hand, the knowledge of its precise value would hardly be of decisive importance 
because the bubble gas parameters are not known precisely. 

Let us describe shortly the proposed scenario. At certain short time interval, r c ~ Ins, when the bubble radius is 
near its minimum, the acceleration of the radius and, therefore, the pressure gradient, assume gigantic values |^] and 
the sign of the pressure gradient is such as to create a strong (depolarizing) electric field directed to the centre of the 
bubble in a thin water layer, ~ l^m, near the bubble surface, because of the flexoelctric effect. What happens then is, 
in effect, screening of this depolarizing field by free electrons. Indeed, during the same time interval the temperature 
of the gas sharply increases up to at least several thousand K, which owing to thermal excitation of the bubble gas and 
water in a thin layer near the bubble surface makes the concentration of free electrons appreciable. The free electrons 
are accelerated by the strong "flexoelectric" field up to energies sufficient to generate additional free electrons as a 
result of the electric breakdown of water. The hot electrons also collide with noble gas atoms dissolved in water and 
transfer them to higher energy excited states. Optical transitions from these states produce light radiation with a 
broad spectrum whose shape is determined mainly by the energy distribution of the hot electrons. The latter has 
the form of a Maxwell distribution with an effective electron temperature which can be very high. This is due to the 
fact that the noble gas atoms have a huge number of excited states with very high probabilities of optical transitions 
between them (see, e.g., pp| , pT| ). Since the radiation occurs in the region of very strong and inhomogeneous electric 
fields the observed radiation spectrum is featureless. The characteristic time of the electric breakdown is much less 
than the characteristic time for the last stage of the collapse r c , i.e. the polarization can continue to change after the 
first breakdown and the electric field can reach the breakdown threshold value more than once. As a result, several 
breakdowns can take place during the time interval r c . During each breakdown the noble gas atoms are excited. An 
important specific feature of the noble gas atoms is the existence of long-living (metastable) excited states with a life 

2l|]). Therefore, once excited the noble gas atom can remain in the 
metastable state for the entire time interval r c ~ Ins, and, possibly, for many periods of the acoustic wave (the period 



time of up to several milliseconds (see, e.g., 1 20 
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is about 30/xs). That means that in the stationary state the longer is the pulse the more is its intensity. 

The paper is organized as follows. In Sec. 2 we discuss the flexoelectric effect in water in more detail than in JTJ] 
and consider the sign of the flexoelectric held, which is crucial for the proposed scenario, at different stages of the 
bubble collapse. The scenario of SBSL is presented in Sec. 3 where the values of relevant parameters necessary for 
realizing the scenario are estimated. In Sec. 4 we discuss the energy distribution functions of hot electrons and the 
breakdown conditions in strong flexoelectric fields as well as the spectrum of SBSL. In Conclusions we summarize the 
results of our theory, discuss their relation to the experimental data, speculate about some possible modifications of 
the scenario and point out some problems to solve. 



II. FLEXOELECTRIC EFFECT AND ELECTRIC FIELD NEAR THE BUBBLE 



We have already mentioned that the proposed mechanism of SBSL is based on the flexoelectric effect, namely the 
appearance of a polarization (electric field) due to gradient of density (pressure) . This effect is not widely known 
and it was discussed in |Q but shortly. Besides, the flexoelectric coeficients are still not determined experimentally 
jl5| . That is why we think it is worthwhile considering in more detail the appearance of electric field due to flexoelectric 
effect . 

The flexoelectric effect is a particular case of a more general phenomenon: appearance of electrical polarization, 
P,as a result of gradient of some scalar quantity, e.g., temperature, concentration of a component, mass density, p. 
This effect should take place in any substance |l6| . As usual, the material equation can be written in several forms. 
In particular, considering polarization as a function of dilatation, u = Ap/ p, and the electric field, E, one has for an 
isotropic medium 

P=aVu+e xE, (1) 

where x is the electric susceptibility and a is one of the flexoelectric coefficients. We shall be interested in the case of 
spherical symmetry and the absence of free charges. In this case D = e E + P = 0. Then from Eqjl] it follows that 

E = -— Vu=- — Vp = /Vp, (2) 

esq eeo 

since u = f3p, where p is the (excess) pressure, /3 is the compressibility. 

To estimate the flexoelectric coefficients it is convenient to consider separately two main mechanisms of polarization: 
charge displacement and dipole ordering. We shall see that in both cases the coefficient of proportionality between 
polarization and gradient of pressure, /, has the same characteristic value. 

The charge displacement polarization is realized, e.g., in ionic crystals. When a gradient of dilatation takes place 
there are less and more compressed regions in each unit cell. The ions of larger radius tend to displace to the less 
compressed region while ions of less radius displace into the opposite direction. Since cations and anions usually 
have different radii their displacements produce polarization. The coefficient a can be estimated as follows K& . 
The maximum ("atomic") gradient of dilatation is equal to l/d, where d is the interatomic distance. This gradient 
has to produce "atomic" polarization P at ~ e/d 2 , where e is the electron charge, when the electric field is absent 
(compensated by some charges), i.e. Ui on ~ e/d . Then from Eq.^| it follows that 

f ^ aion P ^ flL (3) 
seq ds® 

where we have taken into account that e ~ 1 for ionic crystals. 

We shall argue that the same estimation is aplicable for materials with dipole ordering. The dipole ordering 
polarization in the absence of an external electric field but under a gradient of pressure arises because of geometric 
asymmetry of the molecular dipoles. For example, the geometric shape of water molecule is highly asymmetric: the 
negatively charged end is much more compact than the positive one: the negative charge is concentrated in the oxygen 
ion (O ) while the positive charge is shared by two hydrogen ions (H + ) located fairly far from each other. Under 
a pressure gradient a water molecule tries to orient itself in such a way that the oxygen ion would be located in 
the region of higher density while the hydrogens would be located in the region of lesser density. In other words, a 
gradient of pressure leads to orientation of the water molecular dipoles. It is well known that polarization due to 
dipole ordering is much more effective than polarization due to charge displacements. This is reflected in the fact that 
dielectric constants of dipolar materials are, normally, much larger than of non-polar materials. For the same reason 
it is natural to estimate adi P ~ eoti on . Then from Eq.g it follows that for water 
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/ _ adi pP „ fl- „ (4) 

eeo deo N 

where we have taken into accout that for water /3ks 5 ■ 10 _10i ^, d ~ 10~ 10 m. 

For what follows the sign of the flexoelectric coefficients is of importance. As it has been mentioned above the 
negative tip of the water molecule, the oxygen ion, tends to be located in the region of higher pressure. This means 
that the polarization vector is directed opposite to the gradient of pressure, i.e., the flexoelectric coefficient a is 
negative and the coefficient / is positive. 

Let us mention that we rather underestimated the flexoelectric coefficient of water than overestimated. Indeed, 
what we call ai on is replaced, in fact, by a a t defining an "atomic" flexoelectric coefficient that is of the same order 
of magnitude for all substances, i.e. it does not take into account specific features of the substance in question. It is 
natural to expect that a substance consisting of molecules whose electric asymmetry (existence of a dipolar moment) 
is accompanied by a pronounced geometric asymmetry (as it is for water) will exhibit a stronger flexoelectric effect 
than that estimated above. However, while this coefficient is not measured (reliable calculations seem much more 
problematic than measurements) we shall assume the value given by Eq.|| On the contrary, the conclusion reached 
about the sign of the coefficient /, which is crucial for the mechanism discussed in this paper, seems much more 
definite. 

Above we have neglected the conductivity of water. This may seem questionable because the Debye radius of the 
electric field screening in water is comparable with the characteristic size of the strong field region (~ lpm). However, 
this neglect ion is justified because, as a rule, the dielectric relaxation time, td, is much greater than r c ~ Ins, which is 
the maximal characteristic time of the polarization change in our case. Increasing the ionic conductivity of water by 
adding, for example, NaCl, one can, according to our estimations, decrease td down to 0.1ns. In such an electrolyte 
the mechanism of SBSL discussed in this paper might be less effective. 



III. THE SCENARIO: ESTIMATION OF THE MAIN PARAMETERS 

The dynamics of bubble cavitation have been studied in many papers Here we will only discuss the short 
time interval when the bubble radius R is close to its minimum value R c (Fig. la). Fig. la reflects the most essential 
features of Fig. 4 of Ref. || and Fig. 12 of Ref. 0] where experimental results were presented. Within this interval the 
velocity of the bubble surface v — dR/dt reaches its maximum and reaches zero at the point R = R c (Fig. lb) and 
the acceleration a — d 2 R/dt 2 reverses its sign and can achieve huge values (Fig.lc). During the negative acceleration 
period (t ~ ti, see Fig. lc) the gradient of pressure Vp = —pa is directed from the bubble center to its periphery, thus 
according to Eq.|] the flexoelectric depolarizing field has the same direction (Fig. Id). That means that free electrons 
that could be generated in the gas or in water near the bubble surface cannot be accelerated. During the positive 
acceleration period [ti < t < t%, see Fig. le) the situation is opposite: the arising flexoelectric field is directed to 
the bubble center and the generated electrons can be accelerated by the field and produce a breakdown (Fig. le). 
According to experimental data for a bright SBSL the characteristic value of the acceleration is about (10 n -f- 
10 13 ) m/s 2 during an interval, r c = t% — t\ ~ Ins , so the pressure gradient Vp = —pa can reach (10 14 -^10 16 )N/m 3 . 
Note that the gradient of pressure decays with distance from the bubble surface . For purpose of estimations one can 
consider water as incompressible liquid where Vp = (Vp s )i? 2 /r 2 where Vp s is the gradient of pressure at the bubble 
surface. The same estimation for the gradient of pressure can be obtained if we take into account that the pressure is 
about (10 6 -^10 8 ) N/m 2 at the final stage of the bubble collapse and the extension of the high pressure region is about 

It follows from Eqs.[| and J| that for Vp=(10 14 -^10 16 )N/m 3 the value of the flexoelectric field, E ! can reach (10 7 -f- 
10 9 )V/m. According to Eq.|[the electric field decays just as the gradient of pressure, i.e., the flexoelectric field in water 
is given by the formula E ~ E s R 2 /r 2 where E s is the field at the bubble surface, i.e. the radial extension of the strong 
field region is about lpm. The field E can essentially exceed the threshold field of dynamic electric breakdown of 
water E t h which is about (l-^3)-10 8 V/m Q. It does not mean, however, that electric breakdown will occur: a "lucky 
electron" which capable of provoking an avalanche is needed. A similar situation takes place with the laser breakdown 
p8[ . At the same time the conduction electron concentration, n, in water at room temperature is practically zero 
(n ~ 10 -90 cm -3 ): water could be considered as a wide gap amorphous semiconductor with E g ~ 6.5eV fl^ , |T7| , P4 -p4[ . 

The breakdown starts only the moment when the strong flexoelectric field is directed to the bubble center (Fig. If) 
and conduction electrons appear because of the sharp increase of temperature near the bubble surface . The latter 
takes place when the bubble radius R ~ R s which is close to R c ~(0.5-^l)p. At this moment the flexoelectric field 
Ef may be much more greater than E t h so the coefficient of avalanche multiplication of an electron can be practically 
infinite and in this case just several electrons are enough to provoke the breakdown and to screen the field. 
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FIG. 1. Illustration of dynamics of the pulsating bubble and of formation of the flexoelectric field and the electric breakdown 
at the bubble collapse. Time dependence of the bubble radius (a), velocity (b) and acceleration (c) of the bubble surface. The 
flexoelectric polarization at different moments: (d) - at the beginning of the collapse (i < t\), (e) - not very close to the turning 
point, (f) - in the vicinity of the turning point (R — R c ). 

Let us show that even relatively moderate temperatures near the bubble surface (which seem to be generally 
acceptable for the regimes without the shock waves) are quite sufficient to provide a sufficient number of conduction 
electrons to ensure the breakdown. As in Ref. |6| we assume that the gas has the temperature of 7000K. That means 
that the water layer with a thickness of the thermal penetration length, St, has a temperature of about 3000K. For 
5t one has: 
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St = (-t) 1/2 (5) 
pc 

where k, the thermal conductivity of water is about 0.4 J/m s K; c, the specific heat, is about 4-10 3 J/m 3 K; p, the 
density, is 10 3 kg/m 3 . Bearing in mind that the characteristic time for the last stage of the collapse is r c ~lns ||, one 
finds that the thermal penetration length is St ~ 0.1pm. Following |l7],[l8],|4| we will consider water as an amorphous 
semiconductor with a band gap E g ~ 6.5eV and an effective density of states N* ~ 10 21 for T ~ 3000K. The 
equilibrium concentration of conduction electrons in water at T > 3000K is n = N* exp(—E g /2kT) > 3 • 10 15 cm~ 3 . 
Taking into account that the volume of the hot surface layer is approximately V ~ 10~ 12 -j- 10~ 13 cm 3 we find that 
the equilibrium number of conduction electrons in the layer is N > 300 -j- 3000. It is important also to find the time 
needed for the equilibrium concentration of electrons to be established. The thermal ionization frequency (see, e.g., 
{§) is 

v T = (N*av T ) exp(-E g /2kT), (6) 

where a is the cross-section of the free carrier recombination, vt is the thermal velocity of conduction electrons in 
water. Taking into account that a reasonable value of a is about 10~ 15 cm~ 2 and vt ~ 10 7 cm/s we find from Eq.^ 
that for T > 3000K the transient time to the equilibrium tt — v^ 1 < 10~ 8 s. From that follows that during the 
characteristic time for the last stage of the bubble collapse r c ~lns the number of conduction electrons in the layer, 
N, exceeds 30 -j- 300 electrons, which is quite enough to provoke breakdown. 

Note that a similar number of free electrons can be generated by thermal ionization of the gas. Indeed, the assumed 
gas temperature (7000K) is about two times larger and the ionization energy both of the water vapor and Ar, Kr, 
Xe is about 12-^16eV, i.e. also about two times larger than E g . Authors of many works (see ^| and the references 
therein) state that much higher temperatures can be reached because of shock waves forming in the bubble gas during 
the collapse. In principle, one can imagine a situation where high temperatures are not reached (and free electrons 
are not generated) before the shock waves are formed. In this case the shock wave, when it explodes, will start the 
electric breakdown. 

As a result of breakdown the depolarizing flexoelectric field becomes screened. The total transmitted charge in the 
process of the screening is 

Q t = PS = 4ttR 2 s £ E s (7) 

where we have taken into account that in our case D = P + EqE = 0, i.e., P = —EqE. Assuming R s ~ lpm and 
E s = (10 8 -j- 10 9 ) V/m we find that the maximum total number of transmitted electrons is N t = Qt/e ~ (10 5 -j- 10 6 ). 

Note that for considered electric fields the time of the screening (breakdown) , r s , is determined by the time it takes 
for the conduction electrons to cross the region of the strong electric field whose size is Z ~ R s . Thus r s = l/vj, 
where Vd is the electron drift velocity which in strong electric fields saturates at some value Vd ~ 10 5 m/s (see Sec.4) 
Assuming I ~ lpm we find t s ~ lOps. Note that the value of r s is less than the observed SBSL pulse width 

The whole period of positive acceleration, r c , (see Sec. 3) is about two orders of magnitude larger than Tb- T : 
after the breakdown is finished and the depolarization field is screened the polarization continues to change, because 
of change in acceleration. The electric field arises once more and can exceed the breakdown threshold. As a result 
a new breakdown will take place. Such a situation can be repeated several times. Effectively, it manifests itself in 
an increase of the pulse duration which may achieve a fraction of a nanosecond. Therefore, within our scenario the 
greater is the pulse duration the greater is its energy. Analogous interrelation is observed in experiment. |BE6LE7|. 




IV. ENERGY DISTRIBUTION OF HOT ELECTRONS IN THE ELECTRIC FIELD AND SPECTRUM OF 

THE SONOLUMINESCENCE 

We have already mentioned that the electric fields under consideration are very strong, inhomogeneous and change 
over time. However, the characteristic space and time scales of the field change are correspondingly I ~ 1pm and 
t c ~ Ins and they are much more than the relaxation l ength A and relaxation time r £ of the hot electron energy 
which are respectively about (10 -t- 100) A and 10~ 13 s. JT^,|l^Q. Therefore one can consider the local electron energy 
distribution function using well known results for homogeneous static electric fields |^8| Q . 

The considered electric fields are superstrong, i.e. the following condition is valid: 

qEX > hujph = £ph, (8) 

where e p h is the characteristic energy of local oscillations in water which is practically equal to the energy of optical 
phonons in ice, e p h — (80 -j- 100)meV [[17 



18 24]. Using A ~ (10 -j- 100) A one sees that the condition H is satisfied for 
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fields E > 10 7 V/m. The condition | means that in these helds an electron acquires, on average, an energy qEX which 
is much more than e p h- Since in the process of the acceleration an electron generates many phonons the electron 
energy distribution is nearly isotropic in the momentum space. Specifically, in this case the energy distribution of 
hot conduction electrons is approximated with high accuracy, up to energies of electrical breakdown e = E g , by the 
Maxwell function ^fQ: 

/(e^exp^JL) (9 ) 

with an effective electron temperature T e which is determined by the balance equations for energy and momentum of 
the hot electrons 

de e ph v T , . 

_ = eEv d - — = 0, (10) 

— = eE - m e v d — = 0. (11) 

where vt is the effective thermal velocity of electrons. Since for the Maxwell distribution 

3,„ 1 



2 kT e = ^m e vl, (12) 



from Eqsjl^, and it follows that 



Vd 



yje p h/m e (13) 



and 



kT e = (eEX) 2 /3e ph . (14) 



From Eqs [l2] and as well from Eq. || one sees that vx/vdr — \J kT e /e p h = qEX/e p h >> 1. This is precisely the 
condition of validity for the thermalization of the hot electrons and use of the Maxwell distribution. 

Conduction electrons in water form polarons which are conventionally called hydrated electrons ]2^]24|| . However, 
in the strong electric field the polaron states decay and the current carriers are the usual conduction electrons. For 
purpose of estimation we assume that their effective mass m e is close to that of free electrons mo. Then from Eq.[l3] it 
follows that Vd — 10 5 m/s for e p h — lOOmeV. It has already been mentioned that the electric fields under consideration 
are E ~ (10 s -r 10 9 )V/m. It follows from Eq.|| that for such fields kT e ~(H 10)eV, i.e., T e ~ (10 4 10 5 )K for 
A ~ (10 -r 100)A . One sees that the electron temperatures can be several orders higher than the gas temperature in 
the bubble. Recall that for our scenario (see Sec. 3) it sufficient that the gas temperature be about 7000K. 

Now we will discuss the role of noble gases dissolved in water. An important feature of an noble gas atom is the 
existence of metastable states. The life time of the metastable states can reach several milliseconds when the noble 
gas atoms are impurities in solids (see, e.g., [ p0"|]2l| l ) . Since the nearest order in water is essentially the same as in 
solids it would be natural to expect that the life-time of the metastable states for the noble gas atoms in water is not 
less than Ins . Hot electrons not only generate new conduction electrons (the breakdown avalanche) but also excite 
the noble gas atoms into metastable states. Below we will estimate from the experimental data that the characteristic 
value of kT e ~ (2 4- 5)eV. However, a considerable number of electrons may have energies about lOeV and can excite 
the noble gas atoms. 

Note that the energies necessary for excitation to the metastable states, e m , for Xe, Kr, Ar, Ne, He are about 
lOeV and increase monotonically from Xe to He. In other words, it is much easier to excite Xe, Kr and Ar than Ne 
and He. This could be the reason for the increase in SBSL intensity in the series He-Xe, which has been observed 
experimentally |3j3l|. 

As we have mentioned above, the life-time of the metastable states of noble gas atoms in water is expected to 
be fairly long. Thus, once excited a noble gas atom can remain in a metastable state for the entire time interval of 
positive acceleration and multiple breakdowns, r c ~ ins. The major part of electrons has the energies e < e m E g ;they 
collide with these metastable atoms and transfer them to higher excited states. The radiation transitions between 
high-energy excited states of noble gas atoms govern the SBSL spectrum. 

Note that the life-time of the metastable states can exceed the period of the acoustic wave. In this case there will 
be an effect of accumulation of the noble gas atoms in metastable states resulting in the gradual build up of the SBSL 
power during several periods of the acoustic wave. 
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A specific feature of a noble gas is an abundance of excited states with energies higher than those of metastable 
states and with high probabilities of the radiation transitions between them. That is why the optical spectra of the 
noble gas atoms contain many lines Because of the Stark effect in the strong electric field these lines are 

split. Besides, in the active region the strong electric field changes at least several times and the observed radiation 
spectrum is a superposition of the spectra of the atoms in different strong electric fields. In other words one can 
consider the density of the atomic excited states as a constant. 

The probability of a hot electron having an energy e < E g is given by Eq.^ and at every collision it transfers the 
noble gas atom in a metastable state to a state with energy e, the reference point of energy being the energy of the 
metastable state (for our estimations we consider only one metastable state). The concentration of atoms excited 
during the time of a single electric breakdown, i.e. the screening time, r s , to energies within an interval de reads 

£ ds 

dn** = n* n (a ex v T n)T s exp(- — ) • — (15) 

where n* is the number of the noble gas atoms in the metastable state, a ex is the cross section of the impact excitation 
of an atom from the metastable state to a state with energy e. Atoms excited to states with the energy e generally 
do not go to the ground state directly but through intermediate excited states. It reasonable to assume that with a 
high probability they radiate phonons with energy hv ~ e. So the spectral density of the SBSL radiation energy per 
pulse for unit volume can be written as 

~ t hv d(hv) 

P(hv)d(hv) = hvw r n n [a ex VTn)T s exp(--^ r ) — (16) 

where w r is the probability of the spontaneous radiation transition. Taking into account that w r = ^^l^ D 2 where 
D is modulus of matrix element of the dipole moment of the transition p2j we find 

win u/, x 4(2tt) 4 L> 2 , 4 . hv ,d{hv) 

P(hv)d(hv) = n* n {a ex v T n)T s v 4 exp(- — )-^. (17) 

Recall that T e in Eq. [Ij] is a function of the coordenates because T e ~ E 2 (see Eq.[l4| ). Putting, as before, 
E = fWp ~ E s R 2 /r 2 and integrating approximately Eq.|l^ we find the spectral power of SBSL for a single breakdown 
is 

f ~, , 4(2tt) 4 L> 2 , , , hv s 

P = I P(hv,r)dV = 3c3fe <(a e ^ T iV t )T^ 3 exp(- — ) (18) 

where kT es = (eE s X) 2 /3e p h and N t is, practically, the total number of electrons participating in the breakdown. Note 
that the approximative Eq.|l8| is not sensitive to the form of decay of the field, it is important only that E decays 
more steeply than i.e. it is not important that for Vp we use an expression for an incompressible liquid. 

The observed spectra are cut off in the shortwave region because of the absorption of water. This can be taken into 
account by multiplying Eq.|l8| by exp(-a(hv)L) where L is the size of the acoustic resonator (L ~ 2.5cm [|J). Because 
of the Urbach absorption tails fl7|,|23|] the radiation maximum is located considerably below than hv ~ E g ~ 6.5eV\ 



It should be emphasized that the spectrum given by Eq.18 resembles the black-body spectrum but temperature 
is given here by the effective temperature T es of hot electrons near the bubble surface. The value of T es may 
be much higher than the bubble gas temperature. Experimentally observed spectra can be fitted, in the wavelength 
interval (200-^700) nm, to the black-body ones with temperatures (2-^5)10 4 K pj33|. According to EqJlJ such electron 
temperatures are reached at electric fields E s ~ (2 -f- 10)10 8 V/m for values of e p h — O.leV and A ~ (10 100) A 
characteristic for water Jl^j . As follows from EqsJJ and Vp = —pa such fields are realized at accelerations a ~ 
(10 13 -j- 10 14 ) m/s 2 . Similar accelerations are reported in ^j^- 

Integrating Eq.^ multiplied by exp(-a(/iu)L) we find an estimate of the radiation energy for a single breakdown 

W r = Pw~ 1 exp(-a(hv)L)d(hv) ~ (a ex v T N t )T s n* n hv ~ —{a ex R s n*)N t hv (19) 
J v d 

where hv is the characteristic photon energy which is close to the energy of the maximum of the observed spectrum, 
which is about (5 ~ Q)eV in the case of strong electric field of our interest. Recall that in Eq.[l9] r s ~ R s /vd is the 
screening (breakdown) time and Nt — ^i:R 2 s £oE s / e is the number of transmitted electrons in the screening process 
(Sec. 3) and n* is the concentration of noble gas atoms in a metastable state near the bubble surface. 

As we have mentioned before it was shown in ||]4| that because of dissolved gas diffusion and chemical reactions 
an accumulation of noble gas in the bubble takes place while the bubble pulsates and in stationary conditions the 
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gas in the bubble consists almost entirely of noble gas. Therefore one can assume that the concentration of the noble 
gas atoms, n„,near the bubble surface is close to its saturation value. Note that this value has to be estimated for 
atmospheric pressure and ambient temperature rather than for the high pressures and temperatures existing at the 
last stage of the collapse over a very short time. This is why we assume that n n ~ (10 18 — 10 19 ) cm~ 3 . Now we will 
argue that the value of n* can be only one order of magnitude less than n n . Indeed, for kT e ~ (3 -f- 5) eV obtained 
above by fitting the SBSL spectrum to the black-body spectrum, the number of electrons having an energy higher 
than lOeV, Nh may be about 10~ 3 -j- 1CT 2 of the total electron number N t ~ 10 6 i.e., Nh is about 10 3 ~ 10 4 . Such a 
"superhot" electron excites a noble gas atom to the metastable state over the time r ex = (cr ea; UTn„) _1 ~ 10 _12 s for 
a ex ~ 10 _15 cm 2 , vt — 10 6 m/s and n n ~(10 18 -j- 10 19 ) cm -3 . Therefore a "superhot" electron during its participation 
in screening excites T s /r ex noble gas atoms and the total number of the excited atoms is NhT s /r ex ~ 10 4 -j- 10 5 if one 
takes into account that t s ~ I0r ex ~ 10 _11 s. The total number of the noble gas atoms in the region of the strong 
field (~l/im 3 ) is about 10 6 -j- 10 7 , i.e. the percentage of the noble gas atoms excited to the metastable state during 
a single breakdown can be about 10%. Since the life time of the metastable state is much more than r c and during 
this time interval several breakdowns can take place (Sec. 3), one may expect that the number of the noble gas atoms 
excited to the metastable state during the interval r c is only one order of magnitude less than the total number of 
these atoms i.e. n* ~ 10 18 cm~ 3 . If the life time of the metastable states exceeds the acoustic period (~ 30/zs) an 
accumulation of noble gas atoms in the metastable states can also occur during several acoustic periods. 

Now we can estimate the total radiation energy. We have mentioned above that the maximum value of Nt — 10 6 , 
R s ~ 1/im, a ex ~ 10 _15 cm 2 and vt » v<i- Taking into account that the positive acceleration period t c ~ Ins and 
the breakdown time t s is about 0.01ns one can assume that the number of breakdowns is not less than 10. Thus one 
finds from Eq.[l9| that the maximum photon number in the SBSL pulse can be more than 10 7 . This value corresponds 
to the maximum photon number observed experimentally [^],[53 34 1. Note that the radiation energy is a small part of 
the total energy of the flexoelectric field (see Ref |55|). 

V. CONCLUSIONS 

Let us emphasize once more that within the mechanism of SBSL considered the main processes occur in water near 
the bubble surface unlike the most widely discussed mechanisms of SBSL B. Within the framework of this mechanism 
many experimental data about the SBSL can be explained quite naturally. Of course, since the parameters of the 
collapsing bubble are not reliably known when the bubble radius is close to its minimum value we can present no 
more than order-of- magnitude estimations. 

Let us summarize the main results. 

(i) . The minimum duration of the SBSL flash is determined by the single breakdown (screening) time, r s ~ lOps. 
A larger time is possible because of the multiplicity of breakdowns. The maximum duration is limited by t c ~ Ins. 
The longer is the duration the greater is the energy of the flash. This is in agreement with experiments where it was 
found that the pulse duration changes from 30ps to 400ps when its energy increases p|, p6| , p7| | . 

(ii) . According to our estimations the maximum energy of the flash corresponds to 10' -j- 10 8 photons with energy 
(5 -i- 6)eV, which also agrees with the experiment [g 33 3J 



(iii) . Within our scenario the noble gas atoms play an important role. They do not reduce so much the breakdown 
threshold as they govern the radiation process. Abundance of optical transitions in these atoms and inhomogeneous 
broadening because of the Stark effect explain the practically continuous character of the SBSL spectrum. 

(iv) . In agreement with the experiment HQ the theoretical spectrum of the SBSL resembles the black-body 
spectrum but the temperature is given here by the effective temperature of the hot electrons which can be about 
several eV, what corresponds to the observed apparent radiation temperature [§,^3). At the same time the gas 
temperature is not directly related to the radiation temperature and can be considerably less than leV. 

(v) . The increase of influence on SBSL intensity in the noble gas series He-Xe observed in the experiment M is 
connected with the decrease, in this series, of the energy of the lowest metastable state. 

(vi) . The pulse width does not depend on the spectral range of the radiated photons. This also agrees with 
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(vii). The effect of synchronization of the light pulses observed experimentally [[2|jl9|] is explained. 

Note also that within our scenario of SBSL the influence of magnetic fields on SBSL is fairly weak. Their influence 
becomes appreciable when they are high enough to hamper the heating of electrons |56| . This might be the reason 
for the decrease of the SBSL intensity in strong magnetic fields which has been observed experimentally J37|. 

It seems that the considered mechanism of SBSL is especially effective for water because of a lucky coincidence of 
several conditions: 







(i) the strong geometric asymmetry of the water molecule, which accompanies its electric asymmetry, provides 
the needed sign of the flexoelectric coefficient . This sign is such that the electric field has the "correct direction" 
(accelerates electron) over the same (very short) period in which electrons are generated because of the sharp increase 
of the bubble gas temperature. The temporal coincidence of these two cicumstances is a possible reason of the effect 
of synchronization of the light pulses observed experimentally |^,[p| . 

(ii) water is a semiconductor with a relatively narrow band gap (E g = 6.5eV) and sufficienly wide conduction band, 
what is necessary for the heating of the electrons in strong electric field. 

(iii) high solubility of noble gases in water what makes possible high concentrations of noble gas in water near the 
bubble filled by the noble gas accumulated there in the process of the bubble pulsations. 

Fulfillment of these conditions gives the answer to the question why water is the friendliest liquid for SBSL 
Some experimentally found characteristics of SBSL have not been explained in this paper, in particular the de- 
pendence of the energy of the pulse upon the partial pressure of noble gases and water temperature These 
dependences may be not necessarily a manifestation of the radiation processes studied in this paper but they could 
be a consequence of a set of various factors. They may be determined by the hydrodynamics of the pulsating bubble, 
the solubility of noble gases in water and also by the temperature dependence of the flexoelectric coefficient of water 
which still remains, unfortunately, unmeasured. 

Although the main features of SBSL seem to be explainable within our scenario, these explanations are qualitative 
rather than quantitative. We are still a long way from a quantitative theory at present. One of the main aims of the 
paper is to stimulate some experiments that could either support or discard the proposed mechanis of SBSL. 

(i) Measurements of flexoelectric coeficients of water. 

(ii) A detailed study of radiation spectra at electric breakdown of water and ice with different concentrations of 
noble gases. 

(iii) Study of influence of water conductivity on the SBSL intensity. 

After performing these experiments it makes sense to develop the theory further. In particular, it will be necessary 
to study in more detail the kinetics of the screening process and of the excitation of noble gas atoms while taking 
into account the time variation of the distribution function of the hot electrons. This problem should be considered 
together with that of the determination of the spatial distribution of excited noble gas atoms near the bubble surface. 
Diffusion and relaxation of the excited atoms should be considered together, of course, with the processes that lead 
to noble gas accumulation in the bubble |||§]. 

In our estimations we assumed that flexoelectric coefficient of water has its "natural" value. This was sufficient for 
the theoretical estimations to be in agreement with the experimental data. In fact, because of the above-mentioned 
strong geometric asymmetry of the water molecule, the coefficient could be even larger. If this is the case some new 
phenomena will take place including influence of the flexoelectric effect on hydrodynamics of the pulsating bubble 
and X-ray radiation induced by electrons of very high energy. 
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